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A method is presented for predicting 
the total response, in both frequency 
and time, of an unknown linear system 
when only the measured continuous 
wave (cw) magnitude is available. The 
approach is based on approximating the 
square of the measured magnitude by a 
rational function, from which various 
system transfer functions in terms of 
complex frequency are deduced. These 
transfer functions may or may not be at 
minimum phase. The corresponding 
impulse response is then obtained by 
taking the inverse Laplace transform of 
the transfer function. The impulse 
response of the minimum-phase case 
rises faster initially to its first maximum 
than the nonminimum-phase counter- 



parts. This result confirms that, for the 
same cw magnitude response, the accu- 
mulative energy contained in the 
impulse response is the greatest when 
the transfer function is at minimum 
phase. Physical meaning of the energy 
content is also discussed. 
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1. Introduction 

The total response in frequency and time of a 
system to an assumed excitation, whether it is 
continuous wave (cw) or pulsed, is usually unpre- 
dictable if the system involves nonlinear elements. 
Even for linear but complex systems, the task of 
obtaining the total response when the system is 
excited by a source is still formidable because the 
complete system transfer function (amplitude and 
phase) may be unknown. The transfer function is 
defined as the ratio of system output to input in the 
frequency domain. The output and input can be 
voltage, current, electric field, magnetic field, or 
combinations of them. The time response of such a 
system required for assessing its vulnerability to an 
unfriendly electromagnetic environment can be 



determined only by sophisticated time-domain 
measurements or by derivations from the 
frequency-domain amplitude and phase measure- 
ments. Unfortunately, such time-domain measure- 
ments or frequency-domain phase measurements 
require expensive equipment and special consider- 
ations on radiation hazards, regulatory compliance, 
and environmental pollution (if performed out- 
doors). On the other hand, measuring magnitude 
response data of an unknown, complex, linear 
system to cw excitations at low levels, indoors or 
outdoors, is relatively straightforward, less costly, 
and free from compliance and pollution problems. 
Further, if such measured cw magnitude data can 
be processed to deduce a system transfer function, 
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the phase characteristics and time response of such 
a system to a general excitation are then derivable. 

In this paper, we present a method to deduce the 
total response of an unknown, complex, linear 
system from a given set of cw magnitude responses 
only. This is accomplished by approximating the 
square of the measured magnitude curve by a sum 
of ratios of two polynomials with real coefficients. 
Each of these ratios represents a second-order 
rational transfer function for a system with time- 
invariant lumped-constant elements. The exact 
number in the sum is determined by the number of 
resonant frequencies displayed in the measured 
magnitude data. Once the approximation work is 
done, the associated system transfer functions can 
then be obtained by using knowledge available from 
classical network theory. The transfer functions so 
obtained may or may not be at minimum phase. The 
corresponding phase characteristics and impulse 
responses are determined in a straightforward man- 
ner. For the same cw magnitude response, the accu- 
mulative energy associated with the impulse re- 
sponse corresponding to a minimum-phase transfer 
function is always greater than that corresponding 
to nonminimum-phase transfer functions. The deri- 
vations and detailed analysis are presented in Sec. 4. 

The theoretical relationship between amplitude 
and phase of a system with a minimum-phase trans- 
fer function is outlined in Sec. 2. The conventional 
numerical approach for determining the phase 
characteristics and the corresponding time response 
of the linear minimum-phase system from the 
measured cw magnitude data, and the accuracy 
involved in this process are reviewed in Sec. 3. 
Examples are given in Sees. 4 and 5 to demonstrate 
the usefulness of the proposed method developed 
in Sec. 4. Energy contents associated with the given 
system are discussed in Sec. 6. 

2. Theoretical Background 

A stable linear system, however simple or com- 
plex, can be characterized by its transfer function 
H(s), which has no poles in the right half of the 
complex frequency .y -plane. That is, H(s) is analytic 
in Re(s)^0, where Re stands for the "real part of 
[1]. We address only stable systems in this paper, 
because otherwise the system is not well designed, 
and therefore is not useful in application. In addi- 
tion, when the system is made of only time-invariant 
and lumped-constant elements, its transfer function 
is then a rational function ofs (a ratio of two poly- 
nomials with real coefficients) with the degree of 
the numerator polynomial lower than the degree of 



the denominator polynomial. When this transfer 
function is evaluated at s=jo), H(jo)) is then a 
complex function of o), consisting of a real part 
R (o)) and an imaginary part X((o), or a magnitude 
\H(jco)\ and a phase 0((o). That is, 

H{jw)=R{(o)+i X{<o) = \H{ioy)\ e-^"'"), (1) 

where the convention of assigning a minus sign to 
the phase function is used. The magnitude function 
\H{j(t))\ may also be expressed in terms of the 
attenuation function a(w): 

|^0"^)|=e-''(->, ln|^(»|= -a(a)). (2) 

When H{s) is analytic as defined above and the 
system under study is causal [/i(f) = when («0], 
as is usually the case in practice, the real and imag- 
inary parts of H(Jo)) are related by the Hilbert 
transform pair [2], 



R(o>) = ll b>X(y)Kco'-y')]dy. (3a) 

~ Jo 



and 



^(^)=-^f my)Kco'-y^]dy. (3b) 
T* Jo 

In other words, the real and imaginary parts of this 
system are not independent. When one part is 
given either analytically or through measurement, 
the other part can be uniquely determined by 
performing one of integrals shown in Eq. (3). The 
complex transfer function H(j(o) is then com- 
pletely obtained, from which the impulse response 
may be derived. In reality, however, Eq. (3) is not 
useful because we cannot just measure the real or 
imaginary part of the system response to a given 
excitation. 

If H(s), in addition to being analytic and causal, 
also has no zeros in the right half of the .y-plane, 
the transfer function is said to be at minimum 
phase, herein denoted by Hm(s). Under this condi- 
tion, the attenuation function a((o) and phase 
function 0((o) are related by another Hilbert 
transform pair [2, 3, 4], 



.(.)=f/; 



[a(y)/(y'-o>')]dy 



= -% \_MH.{Jy)\i{y'-<o')\Ay, (4a) 
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and 

a(a>) = a(0)-^ \ _Je(y)/y(y'-a>')]dy. (4b) 

From Eq. (4b), we see that the attenuation func- 
tion can be detennined completely from a given 
phase function only when a(0) is also known. But, 
for our application, only Eq. (4a) is required 
because we assume that the magnitude (or attenua- 
tion) function is given by measurement. Once 6((o) 
is detennined from Eq. (4a), the entire complex 
Hm(j(o) can be obtained from Eq. (1) because 
\Hm(J<^)\ or "(<»') is already given. The impulse 
response of this minimum-phase system for / 3^ is 
then calculated by the inverse Fourier transform. 



"(')=^ f. 



Hmijw) e'"" dt . 



(5) 



The system's time response to a general excitation 
can be computed by the convolution integral 
hm(t)*e(t), where e(t) represents an excitation, cw 
or pulse, applied to the system input. The success 
of determining the time response from magnitude 
data is based on the assumption that the system's 
transfer function is at minimum phase. The solu- 
tion of impulse response so obtained constitutes 
the only solution. 

In general, however, there may be multiple solu- 
tions, because other possible transfer functions 
with nonminimum phases giving the same magni- 
tude response may exist. One way of obtaining 
them with our proposed method is to be shown in 
Sec. 4. 

3. Conventional Approach 

Since the improper integral in Eq. (4a) is not 
easy to compute, the conventional approach has 
been to apply a transformation of variables known 
as the Wiener-Lee transform to a(w) or \H(ja))\ to 
obtain the necessary 6(a)). When the Wiener-Lee 
transform [2] 



(0 = - tan(5/2) 



(6) 



is applied, the integration interval ( - <» , oo ) for w 
in Eq. (4a) is transformed into (-ir, ir) for 5. The 
original attenuation function a(<u) and phase 
function 6(a)) will be denoted, after trans- 
formation, respectively as A (S) and T(S). Since 

1 



a(a))=- ln\H(j<o)\=- ^ln[R\a))+X\a))] is 



even function of w and 6(a)) = - tan~\X(a))/R(a))} 
is an odd function of w, their respective transforms 
A (S) and 7(5) will be even and odd functions of 5. 
As such, they may be expanded into Fourier cosine 
and sine series, 

^(5)=flo+fliCOs5+a2COs25 + ...+fl„cosn5+... , 

and (7) 

T(S)=bisinS+b2sm28 + ...+b„sinnS+... , 

where the expansion coefficients are determined by 



flo 



=^fj^'^^'=i\?^'^^'' 



(8a) 



If 2 f 

a„=-\ A(d)cosnddS = - A(S)cosnddd 



and 



b„ =i J T(8)smn8dd =| j T(8)smn8d8 . 

(8b) 

When the system under consideration is causal, the 
expansion coefficients in Eq. (8) are simply related 
by [2] 



b„ = -a„ 



(9) 



an 



Thus, when a(a)) or \H(ja))\ is given, A(8) is 
known. Determination of a„ from Eq. (8a) auto- 
matically yields b„ from Eq. (9), which in turn gives 
7(5) and 6(a)) by means of Eq. (6), and hence the 
complex transfer function H(ja)). The impulse 
response is then obtained from Eq. (5). 

The justification for using the Wiener-Lee trans- 
form and the procedures as outlined above seem 
straightforward. The transform succeeds in con- 
verting the original improper integral in Eq. (4a) to 
a proper integral in Eq. (8a). From the application 
point of view, the important question is then: if the 
integral in Eq. (4a) is difficult to compute before 
the Wiener-Lee transform is applied, is it easier to 
compute fl„ in Eq. (8a) after the Wiener-Lee trans- 
form is used? The answer is most likely negative, 
because the integrand in Eq. (8a) involves compli- 
cated transcendental functions. This explains why, 
in practice, numerical computations are required. 
The entire procedure will then involve: (i) conver- 
sion of the measured data of a(a)) or \H(ja))\ into 
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A (5) by Eq, (6), (ii) numerical calculation of a„ 
frora A\S) by Eq. (8a), (iii) construction of a trans- 
fonned phase function T{S) by including only a 
finite number of tenns in the Fourier sine series 
with b„ = —a„, (iv) conversion of T{S) back to 
d{(i}), (v) determination of the complex transfer 
function H(jci)) based on the given a(6>) and the 
newly constructed 6((o), and (vi) numerical compu- 
tation of h(t) by Eq, (5). The numerical solution so 
obtained constitutes the only answer. Other possi- 
ble solutions with nonminimum phases can never 
be found. In addition, each of the above six steps is 
an approximation, thus exerting doubt about the 
accuracy in the final solution [5, 6]. 

Thus, while the Hilbert transform is useful for 
processing measured cw data directly [7], it may 
not always offer advantage, together with the 
Wiener-Lee transform, for deriving the minimum 
phase. 



4. Alternative but Simpler Approach 

Using passive network theory [8, 9], we can 
deduce a rational transfer function H(s) directly 
and exactly from a squared magnitude function 
\H(ja))\^ expressed as a ratio of two polynomials of 
even order in o), where the order of the numerator 
polynomial is at least two degrees lower than that of 
the denominator polynomial. Thus, if an approxi- 
mate squared magnitude function \H(jci))\^ in such 
a form can be obtained from the measured cw 
magnitude data of an unknown, complex, linear 
system to some excitation, the task of deducing a 
rational transfer function, and subsequently, the 
associated phase function and impulse response (in 
time) is then straightforward. We will show later 
that multiple solutions for systems with the same 
\H(jct))\ are possible. The transfer functions so 
deduced may or may not be at minimum phase. In 
this process, we essentially have assumed that the 
original unknown linear system, which may consist 
of distributed elements and other complexities, can 
be approximated by an equivalent passive network 
system with only time-invariant and lumped- 
constanit elements. The approximation is the only 
one involved in the process. The exact order in the 
final approximate \H(j(o)\^ depends on outstanding 
features in the given cw magnitude data. The most 
important feature displaying a strong resonance at 
a particular frequency can be approximated by a 
simple second-order transfer function. 



4.1 Second-Order Transfer Function 

The second-order transfer function may take 
either of the following two forms: 



H2,=A/(s^+as+b), (10a) 



or 



H^(s)=A(5 +c)/(s^+as +b) , (10b) 

where the parameters /I , a,b, and c are all real. In 
addition, we require 



0<a<2^/b 



(11) 



so that the complex poles are in the left half of the 
J -plane. On the other hand, the parameter c in Eq. 
(10b) may be positive, negative, or 0. When c is 
positive, the zero of the transfer function is also in 
the left half-plane (in fact, on the negative real 
axis). In this case, the transfer function is at mini- 
mum phase. When c is negative, the zero is in the 
right half-plane, and the transfer function is at 
nonminimum phase. When c = 0, the zero is at the 
origin, also constituting a nonminimum-phase case, 
and the dc magnitude response at w = is 0. The 
parameter A is used to match the given maximum 
magnitude response at the resonant frequency. 

The outstanding features associated with the 
second-order transfer functions in Eq. (10) are 
examined in the following analysis. 

4.1.1 Second-Order Transfer Function in the 
Form of Eq. (10a) In this case, we have 



H^a (j(^) =A/(b- 0? +jm ). 
The squared magnitude is given by 



(12a) 



\H^(j<o)\'=H^(j(^)H2ai-jco)=H^is)H^i-s)\. 



\s-jw 



ib-(oY+aW w*-(2b-a^)(o^+b 



:, (12b) 



where w is the only variable. 

Setting the derivative of this squared magnitude 
to yields 6>=0 and (o^=(2b -a^)l2. If 6>flV2, 
6> = gives the location of the minimum, while 
6>^=(2Z? -a^)/2 gives the location of the maximum 
representing the location of the resonant frequency, 
herein designated as 



6>g=(2Z>-fl2)/2>0, 



(13) 
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On the other hand, if aV4<b <aV2, |//(0)p will 
be the maximum. In this case, there is no 
resonant frequency. Thus, if a given magnitude 
curve has a resonant frequency at wo (other than 
0), we require 

b>a^/2, (14) 

which is a stronger condition than that in Eq. (11). 
In terms of wo, the squared magnitude can be 
written as 

\H2.U<or= 4 -1 .,, ■ (15) 

(0 —2<oq<o +b 

We then obtain the maximum 



\H^U<oq)\^ = 



h2 ,,4' 

b —(Oq 



(16) 



which is also nonnegative because of Eq. (13), 
where 6 = a)o + 5a^>wo. 
The relative minimum at « = is given by 

\H2a(0)\=A/b. (17) 

When a) -> 00, |//2a(;a))| -> . 

The half-power points may be defined as the 
frequencies at which the magnitude response of a 
linear system decreases to (1/V^) of the peak 
response. The width between these frequencies 
represents a measure of sharpness of the magni- 
tude response near the resonance. These frequen- 
cies are determined by 



or 



i^2«o"«)r = 5i^2«o'«o)r 

(o*-2<oi(o^+b^=2(b^-<oi)., 
which yields 



(O 



^ = (o^±\/b^ -wS. 



(18) 
(19) 

(20) 



If we denote the half-power frequency on the left 
side of the resonant frequency by wi and that on 
the right side by o^, we have 



(02 = (Oo + '\/b^-(Oo , (Oi = (Oo—'\/b^-(Oo, 



(21) 



<U2-<ui=2V^^~<wo , o}2 + (Oi=2(Oo, and 



Clearly, wl and wi are symmetric with respect to 

! 
). 

The quality factor of the system may be defmed 
Q = «o/(w2 - <oi) , (22) 



(ol 



as 



where W2— «i may be called bandwidth of the 
system. 

Mathematically, we need three conditions to 
determine the three parameters^, a, and b. From 
the application point of view, we can express A, a, 
and b in terms of m, (Oi, (02, and \H(j(Oo)\ from 
Eqs. (13), (16), and (21): 



^ =(«!-«?) |//2,(;wo)|/2, 



b^ = (oi+^((ol-(o^y, 



(23) 



and 



a^=2[ VW+(«2-wi)^(«2 + «i)V4-wg]. 



Thus, when there is only one resonant frequency in 
the measured cw magnitude curve for an unknown 
linear system, such as that shown in Fig. 1, the 
special features such as wo, «i, <02, and \H(j(oo)\ 
can be read from it. The required parameters^ , a , 
and b can then be determined from Eq. (23), 
regardless of whether the relation wi + W2 = 2wo is 
satisfied. The square of the given magnitude may 
be approximately represented by Eq. (12b), and 
the unknown linear system may be represented by 
the second-order transfer function Eq. (10a), 

H2a(.s)=A/(s^+as+b)= ,^^^ ^^ , (24) 



(s+a/2y + p 



where 



p^=b-aV4. 



(25) 



Once this is done, we then obtain the associated 
phase function in accordance with the convention 
used in (1): 

fe,(w) =Yi arg[//2a(-;a))///2a(/«)] 



= tan-\a(o/(b-(o^)]. 



(26) 



Since a and b are positive, fe, (w) varies from to tt 
when (o varies from to 00 . 
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Fig. 1. A squared-magnitude response with only one resonant 
frequency. 

The corresponding impulse response can be 
determined by taking the inverse Laplace trans- 
form [10] of //2a(s) in Eq. (24), 



^2.(0= I e-""^ 



s'mpt, t^O. 



(27) 



The time response R(t) of this unknown linear 
system to a general excitation e(t), cw or pulse, is 
then given by the convolution integral [10], 



R(t)=e(t)*h2u(t). 



(28) 



If indeed the given squared magnitude curve 
with only one resonant frequency shown in Fig. 1 is 
represented by Eq. (12b) and the transfer function 
is represented by Eq. (24), the system under 
consideration will be at minimum phase, because 
this transfer function has no zero in the right half 
of the 5 -plane. In general, however, the system may 
also be at nonminimum phase because the actual 
transfer function could be represented by a 
product of //2a (5) and an all-pass function. In this 
case, the sjrstem transfer function with a nonmini- 
mum phase is given by 



H,(s)=H^(s)H^i(s) , 



(29) 



with the simplest (first order) all-pass function 
defined as 

//aii(5) = (-y-«)/(5 + «), (30) 

where a is a real and positive number. 



Since |//aii(y<u)p= 1, we have 

]//„ Uo>W = \H^ (;^)r \HAJo>)\' = \H2a (Jo>W ■ 

(31) 

This implies that the same squared magnitude 
curve as that shown in Fig. 1 can be represented by 
either //2a (s) or Ha(s). 

In Eq. (30), the zero of Hm(s) in the right half of 
the 5 -plane is the mirror image of the pole in the 
left half of the 5 -plane [8]. In general, the all-pass 
function may be of higher order with more zeros in 
the right half of the 5 -plane and the same number 
of mirror-image poles in the left half of the 
5 -plane. These zeros are not necessarily limited to 
the real axis. They can take complex-conjugate 
pairs. If we restrict ourself, for the time being, to 
the first-order all-pass function given in Eq. (30), 
we can make fiirther analysis. The impulse 
response for this nonminimum-phase transfer 
function may be determined from Eq. (29) by 
convolution integral [11], 



= /i2a(0-2a/i2a(0*e"'", 



(32) 



where £"Ms the inverse Laplace transform. 

The impulse response can also be obtained by 
taking partial fractions of Eq. (29). That is. 



//n(s) = 



A(s —a) 



(s^+as+b)(s+a) 



^s +as+b s+ai 



-A 



p(s+a/2)+q-ap/2 p 
(s+a/iy + P^ 

where 

p=2a/(b + a^—aa), and 

q = (b — a^ + aa)/(b + a^ - aa) . 

The impulse response is then 



-^}. (33) 



(34) 



/i„(0=^{e-'"'='I> cosiSf + ^-^sini3f]-/» e-j, 

(35) 

which is the same as Eq. (32) after the convolution 
is performed. 
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Both the minimum-phase impulse response h^ait) 
in Eq. (27) and nonminimum-phase impulse 
response hn(t) in Eq. (35) vanish at / =0. This can 
be confirmed by the initial-value theorem [10]. 
More detailed behaviors of Eqs. (27) and (35) can 
be learned by examining their first time derivatives, 



hL (0 =^ e-'/^E/g cos /a -| sin /»] ., 



(36a) 



and 



hl(t) =A\(g -ap) cos /» 



-[ 



a(2q-ap) 



4)8 



+ pp 



I sin /3/| 



.-a//2 



+Apa e" 



(36b) 



Clearly, att=0,hia (0) =/i„' (0) =A . This means 
that both h2a(t) and hn(t) arise from with the 
same starting rate. However, shortly afterward, 
they increase with different rates. At 
t = £(l>£>0), we may expand those functions 
involved and keep the fi'rst-order term to obtain 



and 



/i2'«(e)«^(l-ae/2)2«^(l-fle) (37a) 



/ii(e)«/l[l-(a+a)e]. (37b) 



Since both a and a are real and positive, we con- 
clude that 



hi(e)<hia(e) . (38a) 

By similar steps we can also show that 

hn(e)<h2a(e). (38b) 

The relations in Eqs. (38a) and (38b) imply that 
shortly after the system is excited by a source, 
h2a(t) associated with the minimum-phase system 
increases with a greater rate than the non- 
minimum-phase counterpart hn(t). Although this 
point is drawn from a special case (first-order 
all-pass function), it can be generalized to higher 
orders even though the algebraic derivations are 



much more involved. Equations (38a) and (38b) 
also mean that more energy accepted by the 
minimum-phase system is concentrated at the 
beginning (/ = + ) of the excitation than the non- 
minimum-phase system [7, 12]. Details on energy 
consideration are found in Sec. 6. This observation 
is very important from the standpoint of electro- 
magnetic interferences (EMI). If the minimum- 
phase system can survive the initial impact due to 
an unwanted external source, a nonminimum- 
phase system can also survive it. The minimum- 
phase system may be considered the worst case as 
far as the initial impact due to an unwanted signal 
is concerned. From the design point of view, if a 
system is minimum phase, the designer may wish to 
convert it to nonminimum phase by adding an all- 
pass network to reduce initial EMI impact. 

From Eq. (36a) we know that the first maximum 
of /i 2a (0 occurs at /mo, which is the smallest root of 



tan /3/ = 2/3/fl 



(39) 



After reaching its first peak at /mo, the impulse 
response /i2<i(0 varies sinusoidally with a decay 
rate of a/2 and with a period of )3. Although it is 
not as straightforward to determine the exact loca- 
tion of the first maximum for hn(t) by setting 
Eq. (36b) to 0, we know that it also varies sinu- 
soidally with the same period ^ but decays with a 
different rate because of the extra term exp(— a/ ). 
A numerical example is here presented to illustrate 
this point. 

Example 1. Suppose that the square of a "mea- 
sured" cw magnitude curve can be represented by 



f(6)) = e-('^-^)' 



(40) 



For this example, the resonant frequency occurs at 
ft)o = 2. The half-power frequencies are <ui =3.1674, 
and 0)2 = 4.8326. The bandwidth is given by 
0)2-0)1 = 0.4185, and g =4.7780. Since 
o)i + o)2 = 2o)o, we can use H2a(s). The required 
parameters can be obtained from Eq. (23) as: 
a =0.4141, b =4.0857, and A =0.8326 (carried to 
4 digits). The approximate squared magnitude is 
then 



W:.u<or = 



0.6931 



0)^-80)2+16.6931 



(41) 



By presenting the numerical results in Table 1, we 
see that the given curve in Eq. (40) and its approx- 
imation in Eq. (41) indeed match at o)i, 0)2, and o)o. 
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Table 1. Approximation of the given function in Eq. (40) by the 
squared-magnitude function in Eq. (41) 



P{->) 



\H^.(J(o)Y 



Error 



0.00 


0.0000 


0.0415 


0.0415 


0.50 


0.0000 


0.0470 


0.0470 


1.00 


0.0001 


0.0715 


0.0714 


1.50 


0.0468 


0.1846 


0.1378 


1.75 


0.4152 


0.4409 


0.0257 


1.80 


0.5612 


0.5454 


-0.0158 


1.85 


0.7164 


0.6752 


-0.0413 


1.90 


0.8589 


0.8200 


-0.0389 


1.95 


0.9617 


0.9467 


-0.0150 


2.00 


1.0000 


1.0000 


0.0000 


2.05 


0.9598 


0.9441 


-0.0157 


2.10 


0.8453 


0.8048 


-0.0405 


Z15 


0.6787 


0.6414 


-0.0373 


Z20 


0.4938 


0.4955 


0.0017 


Z25 


0.3234 


0.3804 


0.0570 


2.50 


0.0063 


0.1204 


0.1141 


3.00 


0.0000 


0.0270 


0.0270 



After obtaining Eq, (41), we have the following 
simplest solution to represent the transfer function 
with a minimum phase: 

//a,(5) = 0.8326/(5^+0.414l5 +4.0857; . (42) 

The associated phase function and impulse 
response are then respectively 

0m(<a) = tan-'[a<u/(6 - <u^)] 

= tan-'[0.4141w/(4.0857-w2)] (43) 
and 

/i2a(0 = 0-4141 e-''-^°'sin(2.0107O, ^>0. (44) 

For the same |//a,(y<u)p obtained in Eq. (41), we 
could also have a transfer function with a non- 
minimum phase by including an all-pass function, 
say, the first order with a = \. We then have 



H^is) = 



0.8326(5 - 1) 



(s^+0A141s +4.0857)is +1) 



(45) 



The corresponding phase function and impulse 
function (after taking partial fractions) are 



0n(<u) = 0m(w) + tan-'(w) + Tr, 
with 0ra(<u) given in Eq. (43) and 



(46) 



h„(t) = [0.3564 cos(2.01070 + 0.2735 
sin(2.01070] e""^""* 
- 0.3564 e-', t^O. 



(47) 



The impulse responses obtained in Eq. (44) and 
Eq. (47) are plotted in Fig. 2 to confirm the con- 
clusions in Eqs. (38a) and (38b). Thus, more 
energy is concentrated in /la, (f ) than in /in (0 at the 
beginning of excitation. 




Fig. 2. Impulse responses of the minimum-phase and nonmini- 
mum-phase systems with their given squared-magnitude cw 
response in Eq. (40) and the approximate squared magnitude in 
Eq. (41). 



The solutions for da(o)) and ha(t) are not unique 
because they depend on the choice of specific 
all-pass functions. 

We note from Eq. (23) that the parameter a 
decreases with the bandwidth (<U2-<ui) or is 
inversely proportional to Q, and that the parame- 
ter b is primarily determined by cdq. These parame- 
ters decide respectively, in turn, the decay rate and 
period of variations of the impulse response. 

4.1.2 Second-Order Transfer Function Taking 
the Form of Eq. (10b) For this case, we have 



H2b(j0l)) = 



and 



\H^(}o>r = 



A(c+joi)) 
b-w^+jcDO 



A\cH<o^) 



(48) 



<o'-(2b-a^)(o^+b^ 



(49) 
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Now there are three solutions for d\H2biJci))\V 
da) = 0. One of these, w = 0, is the minimum with 
\H(0)\'^=A^c^/b\ The other solution gives the res- 
onant frequency. 



m 



= -c^ + \^c*+b^ + 2bc^-a^c\ 



(50) 



which is always greater than (as it should be) in 
view of Eq. (14). A third solution takes the same 
form as Eq. (50) except with a negative sign in 
front of the square-root sign. This third solution is 
obviously nonphysical. 

The maximum value of the squared magnitude at 
Wo is, after substitution of Eq. (50) into Eq. (49) 
and some algebraic simplification, given by 

\Ha(jo>o)\^=AV(2a>i-7b+a^). (51) 

Again, this is the only maximum. The general vari- 
ation of |//2i(;w)|^ is similar to that in Fig. 1. 

The half-power frequencies wi and cd2 are deter- 
mined by 



a)^ + c^ 



a)*-(2b-a^)a)^+b^ 4(a)i-b) + 2a' 



or 



a)U(7b-a^-4(oi)a)^+b^-c\4c^ + 2a^-4b) = 0. 

(52) 
Instead of solving for cai and 0)2, we note that 

a)^i+a)l = 4a)i + a^-2b, (53) 



and 



(o\<al=b''-c\4(ol + 2a''-4b). (54) 



The above derivations represent analysis for the 
particular transfer function H-2b{s). From the 
application viewpoint, we can express a,b, and c in 
terms of wo, (ai, and (a^ by using Eqs. (50), (53), 
and (54), 

b^ = 2<aa — <x)i<x)2, 

c2 = (W-w?w|)/(w? + w|-2W) , (55) 

a'^ = <al-\-o)2-4o)l + '2b. 

Thus, when wo, wi, and <xi2 are read from a mea- 
sured magnitude curve, we can easily determine 
the required parameters b,c, and a from Eq. (55). 
Since a^ b^, and the denominator of c^ [which is 



equal to the denominator of Eq. (51)] are all posi- 
tive, and c ^ itself is nonnegative, we require 



<x)\ + (i)2> 2(1)0 and (00^(01(02. 



(56) 



In addition, the constant factor A can be deter- 
mined from Eq. (51), 

A^ = (2a>i-2b+a^)\HU(Oo)W (57) 

where |//(;wo)|^ can also be obtained from the 
given magnitude curve. 

For the special case c=0 when 0)0 = 0)1(02, we 
have 

b =0)0 = 0)1(1)2, 

a^ = 0)i+0)2 — 2Wo = W 1 -I- W2 - 2Wi W2 = (W2 - ft)i )^ , 

[OT a =o)2-o)i], (58) 

and 

A^=a^\H(jo)o)\^ , or ^ =a |//(M)| . 

Here, the parameter a is controlled solely by the 
system's bandwidth, and b depends only on the 
resonant fi-equency. Also, for this special case, wo is 
the geometrical mean of wi and ^2 . 

One unique feature associated with this case 
[Eq. (10b)] is that once an approximate squared 
magnitude in the form of Eq. (49) is obtained from 
the given magnitude curve by the procedures thus 
outlined, the solutions for the transfer function are 
not unique. One of the obvious solutions, herein 
designated as Hm(s) takes the same expression 
given in Eq. (10b), 



A(s+c) 
s +as +b 



(59) 



where c is taken as a positive number from 
Eq. (55). In this case, the transfer function is in 
minimum phase. 
The other solution is 



A(c-s) 
s^+as+b 



(60) 



which represents a nonminimum-phase transfer 
function. 

Clearly, we have \H„,(jo))\^ = \H„(jo))\^ . 

The associated phase function for Hm(s) can be 
obtained directly from Hmijo)) as 

0ra(w) = tan-'[aW(fe -w^)]-tan-'Cw/c). (61) 
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The second term becomes ttH when c = 0. Since 
the first term varies from to tt and the second 
term varies from to ttH as a changes from to 
00, the range of variation for 6m{<>>) is hence from 
to 7r/2. 
The phase function for Ha{s) is 

0„(6)) = tan"'[fl<u/(6 - <u^)] + tan" 'Cw/c), (62) 

whose range of variation is from to iirll. 

According to the Hilbert transform given in 
Eq. (4a), the phase function 6m{o>) may also be 
obtained from the deduced squared magnitude, 

0".H=-^ f [ln|//™(;>')|/(y^-6,^)]dy 

TV J — 00 

^ r ln(y^+c^)dy 



+ 



I 

- r 



\r)\y*-(2b-a^)y^+b^^dy 



(63) 



Because the integrands are even functions of y, we 
have 



IT J( 



+ 



IT 



0) 

TT 



f 

JO 

f 

Jo 



y'-co' 


ln(y^+m^)dy 


y'-a>' 


In(y2+rt2)dy 


y'-co' ' 



(64) 



where we have broken the last integral in Eq. (63) 
into two parts with 

mHrt2=-2A+fl^ and mV=b'' . (65) 

Since a known definite integral in the form of 
Eq. (64) is available [13], 



Jo 



hV-k^ 



= ^tm-'(gk/fh), f,g,h,k>Q, (66) 



we easily, after comparing Eqs. (64) and (66), iden- 
tifyg=A =l,k = (o,f=c,m, andn respectively for 
the first, second, and third integrals in Eq. (64), 
and thus obtain 

6m{(o)= -tan"'(6)/c) + tan"V'<>/'M) + tan"'('W«), 
= -tan""(6)/c) + tan"'[fl6)/(& -a^) (67) 

which is identical to Eq. (61). The last step in 
Eq. (67) is accomplished by combining tan~'(<u/»i) 
and tan"'(<u/rt) after using the relations in Eq. (65) 
and noting the requirements in Eq. (66). The phase 
function for Ha(J(o), however, cannot be obtained 
from the Hilbert transform. 

The derivation of Eq. (67) from |//ni(;<u)P is 
exact. In general, there is no further approximation 
involved once an approximate squared magnitude 
in the form of a ratio of two polynomials in even 
orders of (o is deduced from a measured cw magni- 
tude curve. In fact, the phase of the minimum- 
phase transfer function can be obtained directly 
from this deduced squared magnitude with the 
help of Eq. (4a) by the method proposed here. The 
type of integral formula given in Eq. (66) together 
with the Hilbert transform Eq. (4a) can also be 
applied to H^(s) in Eq. (24) to obtain the same 
phase function fe (<u) in Eq. (26). 

The corresponding impulse responses are deter- 
mined from 



H„,(s)=A 



and 



H.(s)=A 



(s+c) s+a/2+c-a/2 

s^+as+b (s+a/2y + p^ ^ ' 



c+a/2-(s+a/2) 
(s+a/2y + p^ 



(69) 



yielding respectively 



c —a/2 
hm(t) = A e"""^ [cos pt + sin pt], t^O, 



)8 



(70) 



and 

h„(t) = -A c-'"^ [cos pt - ^ sin ^t], t^O, 

(71) 
where /3 is given in Eq. (25). 

These impulse functions are still, basically, a 
sinusoidal function with a decay rate of a/2, even 
though the form is little more complicated than 
that of /!2a(0 given in Eq. (27). Again, the decay 
rate is related primarily to the bandwidth, and the 
period of variations to the resonant frequency. 
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From Eqs. (70) and (71), we see that hm(0 + ) 
=A [= - hn(0 + )], which is different from the 
previous case where /i2fl(0 + )=0. Both impulse 
responses start (/ = 0) at the same magnitude A . 
Since the parameter c can be greater than a/2, at 
/ =0+ , hm(t) may increase and reach the maximum 
before falling to its first zero. On the other hand, 
\hn(t)\ always deaeases beginning at / =0+ . The first 
zero oihm(t) is determined by the smallest root of 

tanj3/,„= - p/(c-a/2), (72) 

while that of ha(t) is by the smallest root of 

tan )3/„ = p/(c +a/2) . (73) 

A comparison of Eqs. (72) with (73) indicates, 
regardless of the relative values of c and a/2, that 
the first zero oihn(t), called /„o, is no greater than 
that oihm(t), called /mo. That is. 



tnO^ tJO . 



(74) 



The equality sign holds when c = 0. 
Equation (74) implies that t^ is closer to / = than 
is/mo. Thus, the beamwidth of the impulse response 
hm(t) associated with the minimum-phase transfer 
function is wider than the beamwidth of \hn(t)\ 
associated with the nonminimum-phase transfer 
function with the same squared-magnitude func- 
tion. This, in turn, means that more energy is 
concentrated, at the beginning of excitation, in the 
minimum-phase system than in the nonminimum- 
phase system. Once again, the minimum-phase 
system may be considered as the worst case as far as 
the initial impact of the system by the interference 
source is concerned. Another example is presented 
for illustration. 

Example 2. Suppose that the square of a given 
"measured" magnitude can be represented by a 
shifted Gaussian function. 



/='(&)) = 4 e-^t"-^)' 



(75) 



The resonant and two half-power frequencies are: 
ft)o=3, ft)i =2.4113, and ft)2 = 3.5887, the system 
bandwidth is 1.1774, and Q =2.5480. The maximum 
at ft)o is 4. Since wi-l-wi is not equal to 2ft)o in this 
case, we wish to approximate Eq. (75) by |//2t (jo})\^ 
given in Eq. (49). Using Eq. (55) we obtain the 
required parameters: fl^ = 1.3606, &^ = 87.1181, 
0=^ = 8.8267, and yl^' =2.7725. Then 



\H^(M' = 



2.7725(ft)' + 8.8267) 
ft)''-17.3069ft)=' -I- 87.1181 



which is computed together with /^(w) in Table 2 
to show the quality of approximation. 



Table 2. Approximation of tlie given function in Eq. (75) by the 
squared-magnitude function in Eq. (76) 



(0 


PM 


W^U'oW 


Error 


0.00 


0.0000 


0.2809 


0.2809 


0.25 


0.0000 


0.2864 


0.2864 


0.50 


0.0000 


0.3037 


0.3037 


0.75 


0.0002 


0.3350 


0.3348 


1.00 


0.0013 


0.3848 


0.3835 


1.25 


0.0087 


0.4607 


0.4520 


1.50 


0.0444 


0.5768 


0.5324 


1.75 


0.1757 


0.7579 


0.5822 


2.00 


0.5413 


1.0493 


0J080 


2.25 


1.2986 


1.5323 


0.2337 


2.50 


2.4261 


2.3206 


-0.1055 


2.75 


3.5300 


3.3844 


-0.1456 


3.00 


4.0000 


4.0000 


0.0000 


3.25 


3.5300 


3.3851 


-0.1449 


3.50 


2.4261 


2.3215 


-0.1046 


3.75 


1.2986 


1.5294 


0.2308 


4.00 


0.5413 


1.0396 


0.4983 


4.25 


0.1757 


0.7398 


0.5641 


4.50 


0.0444 


0.5495 


0.5051 


4.75 


0.0087 


0.4231 


0.4144 


5.00 


0.0013 


0.3356 


0.3343 


5.25 


0.0002 


0.2728 


0.2726 


5.50 


0.0000 


0.2264 


0.2264 



The computation is carried out only from w = to 
(0 = 5.50. The approximation around the important 
region near w = ft)o=3 is very good, while that near 
the two ends ((o=0 and (o=5.50) is marginal. 
Physically, however, the less accurate results near 
the two frequency ends are of secondary impor- 
tance. 

If we deal with \H:2b(j(o)\^ alone without includ- 
ing extra all-pass functions, we have two possible so- 
lutions for the transfer function, one with a 
minimum phase and the other with a nonminimum 
phase. 



//m(5) = 



and 



(76) 



1.6651(5-1-2.9710) 
5 ='-1-1.16649-1-9.3337 

1.6651(5-1-0.5832-1-2.3878) 
(5 -1-0.5832)=' -1-2.9989=' 

1.6651(2.9710-5) 
5 ='-1-1. 16645 -1-9.3337 

1.6651[3.5542- (5 -1-0.5832)] 
(5 -I- 0.5832)=' -I- 2.9989=' 



(77) 



(78) 
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The corresponding phases are 

0„(6>) = tan-'[1.16646>/ (9.3337- w^')] 

-tan-'('6>/2.9710) (79) 

and 

e„((o) = tan-'[1.16646>/ (9.3337 -6>^)] 

+ tan-V6>/2.9710) (80) 

The impulse responses are 



and 



;i™(0 = 1.6651[cos(2.99890 

+ 0.7962 sin(2.99890] e""-^"^ 

;i„(0= -1.6651[cos(2.99890 

- 1.1852 sin(2.99890] e-'-^^^^ 



(81) 



(82) 



Thus, this method not only yields a solution with 
a minimum-phase transfer function such as by the 
conventional numerical method described in Sec. 3, 
but also gives other possible solutions with non- 
minimum-phase transfer functions. 

Figure 3 shows hm(t) and \ha(t)\ for comparison 
purpose. Evidently, both hm(t) and \hn(t)^ begin at 
/ = with a magnitude of 1.6651. Then hm(t) 
increases to its maximum of 1.9030 at f =0.1602, 
and —ha(t) starts to decrease. The first zero of 
h„(t) is at f = 0.2337, while that of hm(.t) is at 
f =0.7480. Also, |/i„(OUa < |/ini(Olmax in this case. 
More energy is concentrated in hm(t) than in ha(t) 
near t=0. 



\ 












\ 






^ 


\, 




w% 






y^ 






V 













3 
Time, s 



Fig, 3. Impulse responses with the given cw squared-magnitude 
response in Eq. (75) and the approximate squared magnitude in 
Eq.(76). 



Comparing Eqs. (77) and (78), we may also ex- 
press Ha{s) in terms of Hm{s) and a first-order 
all-pass function: 



H,(s)=H4s) 



2.9710 -J 
2.9710+j* 



(83) 



In general, the solution obtained in Eqs. (69) 
and (71) is not the only one with a nonminimum 
phase. In fact, the same given magnitude curve can 
also be represented by the product of Hm(,s) [or 
Ha(s)] and additional all-pass functions. The solu- 
tions depend on the choice of these extra all-pass 
functions. 

4.2 First-Order Transfer Function 

The two types of the second-order transfer func- 
tion analyzed in Sec. 4.1 are the most important 
ones which can be used to approximate a measured 
magnitude curve with only one resonant frequency 
not at 6> = 0. For studies of radiated susceptibility, 
these may be sufficient because it is difficult or 
meaningless for an antenna to measure the inter- 
ference response of linear systems at (o = 0. How- 
ever, when the dc interference is also possible, in 
addition to the cw interferences, for some practical 
systems, the response at (o = may constitute a 
relative maximum. To cover this case, we can 
approximate this part of the given measured 
magnitude curve by a squared magnitude corre- 
sponding to the first-order stable transfer function 



Hi(s)=A/(s+a), 



(84) 



where A and a are real and positive. Its squared 
magnitude is 



\Hi{jo>)\'=AV{o>'+a'). 



(85) 



Obviously, its only maximum occurs at (o = 0, with 
\Hi(0)\^ = (/i/ay . A representative curve for 
\Hi(jco)\^ is shown in Fig. 4. 
The half-power frequency may be determined by 



l/(a>'+a') = y(2a^). 



(86) 



which yields only one solution (02= a (the other 
half-power frequency (oi=-a has no physical 
meaning). The bandwidth in this case is just 2(02. 
The associated phase function can be obtained 
either from Hi{j(o) or from the Hilbert transform, 



0i(w) = tan"'(w/«). 



(87) 



which varies from to Tr/2 as o) varies from to » . 
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Fig. 4. A squared-magnitude response corresponding to the 
first-order transfer function, normalized Eq. (85). 



The corresponding impulse response is 
hi(t)=Ae-', t^O. 



(88) 



The decay rate in this case equals numerically the 
half-power frequency. 

43 More General Case 

For a more general case where the measured 
magnitude curve has N>1 distinct resonant 
frequencies, it may be approximated with our pro- 
posed method by a sum of terms of the type 
\H2a(Joi)\^ or |i/2*(;co)|^ discussed in Sec. 4.1. To 
simplify the notation, let us temporarily drop the 
subscript a and b while maintaining the subscript 2 
to indicate the order of the transfer function being 
considered. The approximated squared magnitude 
will then take the form: 



mj<ow = i i//2,o-co)r 



(89) 



i"l 



with the required parameters a, b,A, and possibly 
c in each |//2(;c<))|^ to be determined by the out- 
standing features associated with each resonant 
frequency. Should there also be a relative maxi- 
mum at 6) = 0, another term in the form of Eq. (85) 
for a first-order transfer function may be added to 
Eq. (89). Once this approximation is accomplished, 
the system transfer functions can then be deduced 
by the classical method [8]. We can then determine 
from these transfer functions the corresponding 
phase functions and impulse responses to give the 



complete characteristics of the unknown linear 
system. Another example is given below to illus- 
trate this point. 

Example 3. Suppose that the square of a mea- 
sured cw magnitude curve can be represented by a 
sum of two mathematical expressions 



f\<o)=mco)+fi(<o). 



where 



/?(6,) = e 



_ ~-2lD 



(90) 



(91) 



is used to simulate a possible maximum at w = 0, 
and 



fU(o) = 4e-^<---^^' 



(92) 



is simulated for a possible resonant frequenqf at 
6)0=3. The expression in Eq. (92) is the same as 
that presented in example 2. 

Since the maximum of/f(c<)) occurs at w = with 
/i(0) = l, and the half-power frequency is 
6)2 = 0.3466, we approximate fi(co) by a linear 
system with lumped-constant elements represented 
by a squared-magnitude function in the form of 
Eq. (85), 



|//,(;6))P = 



0.3466' 



6)2 + 0.3466'' 



(93) 



which is plotted in Fig. 5 together with fi((o) to 
show the approximation involved. 
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Fig. 5. Comparison of functions given in Eqs. (91) and (93). 
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For /l(w), we have m = 3 with /2(3) = 4. 
The half-power frequencies (wi =2.4113 and 
6)2 = 3.5887) have been obtained in example 2. 
Even though the relation a)o = a>ia>2 is not satisfied, 
this time we choose c=0 in H^b^s) to make a 
better approximation at w = 0. The other required 
parameters are determined from Eqs. (55) and 
(57):fe=w§ = 9, fl=' = w?+w|-2wo^ = 0.6931, fl = 
0.8325, and ^4=^= 0^1(3) = 2.7725. We then have 
the following to approximate the given /2(w): 



\H^(J<ot = 



2.7725w^ 



w''-17.3069w=' + 8r 



(94) 



which is plotted in Fig. 6 together with the given 
/l(w) for comparison purpose. The approximation 
shown here may be compared with that shown in 
Table 2, where c is chosen not equal to 0. 
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Fig. 6. Comparison of functions given in Eqs. (92) and (94). 



The total squared magnitude to approximate 
/^(w) in Eq. (90) is then given by the sum of ex- 
pressions in Eqs. (93) and (94): 



0.1201 2.7725w^ 

\H{j<,)Y=—— + 



0) 



^-F0.1201 w''-17.3069w^+81 



2.8926(6)^ - 0.60366) ' + 3.3639) 
" (6)'+0.1201)(6)''-17.30696)'+81) 



= H{s)H{-s)\s.j. , (95) 

which is plotted in Fig. 7 together with p(o)). 
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Fig. 7. Comparison of the functions given in Eqs. (90) and (95). 

From Eq. (95) we can apply the classical method 
in network theory [8] to extract the transfer func- 
tions as follows: 

2.8926(j * + 0.6036iy ^ + 3.3639) 
H(s)H(-s) - (0.1201 -j')(j*+17.3069j' + 81) 



2.8926A^(j) 
D(s) ' 



(96) 



where 



N(s) = (j^ + 1.7506iy + 1.8341) 
(j^-1.7506iy + 1.8341), 
and 



(97) 



D(s) = (0.3466+5')(0.3466-5')x 

(s ^ + 0.8325^ + 9)(5' ^^ - 0.8325^ + 9) . (98) 

From the formats specifically expressed in Eqs. 
(97) and (98), we can assign appropriate factors to 
H(s) and H(-s). Obviously, the factors with 
positive signs in Eq. (98) have to be assigned to the 
denominator of H(s) because the poles are 
required to be in the left half of the s-plane for the 
linear system to be stable. The remaining factors in 
Eq. (98) with negative signs belong to the denomi- 
nator of H(-s). However, either factor in Eq. (97) 
can be assigned to the numerator of H(s) because 
the zeros can be in the left-half or right-half plane. 
If the zeros are in the left half-plane, the system is 
at minimum phase. If they are in the right half- 
plane, the system is at nonminimum phase. Thus, 
we obtain the minimum-phase system, 
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1.7008(^'+ 1.7506? + 1.8341) 
ffm{s) - (^2+o.3466)(s2 + 0.8325s+9) 

= 1.7008[0.1526/(s +0.3466) 

+ (0.8474s + 1 .3299)/(s ^ + 0.8325s + 9)] , 



(99) 
and the nonminimum-phase system, 

1.7008(j^- 1.7506? + 1.8341) 
"^^^ ~ (s+0.3466)(s2+0.8325s+9) 

= 1.7008[0.2900/(s+ 0.3466) 

+ (0.7100$ -2.2381)/(s 2+ 0.8325s +9)] . 

(100) 

We can verify that \Hm(J(o)\^=\Ho(J<o)\'^ and that 
they both equal the squared magnitude in Eq. (95). 
Now, the nonminimum-phase transfer function in 
Eq. (100) can also be expressed in terms of the 
minimum-phase transfer function in Eq. (99) and 
an all-pass function 



H,(s)=H.,is)H,n(s), 



(101) 



where 



H^is) = (s2-i.7506y + 1.8341)/ 
(s2-(-l.7506y + 1.8341), 



(102) 



is a second-order all-pass function with the 
complex-pair zeros in the right half-plane as mirror 
images of the poles in the left half-plane. 
The associated phase functions are respectively 



0m(w)=0i(w)+02(w)-03(w) (103) 



and 



0„(w) = 0,(w) + 02(w) + 03(w), (104) 

where 0i and 62 are due to the denominator factors 
in Eq. (99) and 63 is due to the numerator in 
Eq. (99). They are: 

0i(w) = tan-'(w/0.3466), 

62(0)) = tan- '[0.8325^/(9-0)2)], 



and 



The minimum phase 6m((o) in Eq. (103) can also be 
obtained from Eq. (95) with the Hilbert transform, 
as demonstrated before. 
The impulse responses are: 

/i„(0/l-7008 = 0.1526 e-o-^^ 

+ [0.8474 cos(2.97100 

+ 0.3289 sin(2.9710O e""""*^ (106) 

and 

/i„(O/l-7008 = 0.2900 e-"-^^ 

+ [0.7100 cos(2.971*) 

- 0.8528 sin(2.971*) e-"""*^ (107) 

where both hm(t) and ha(t) are normalized with 
respect to the common constant factor 1.7008. 

The normalized impulse responses in Eqs. (106) 
and (107) are plotted in Fig. 8. Here, they both 
start with 1 at t=0. The normalized ha(t) 
decreases much faster than the normalized hm(t). 
Thus, ha(t) has a narrower beamwidth than hm(t), 
or less energy is concentrated initially with ha(t) 
than with /im(0- 
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03(w) = tan-'[1.7506w/(1.8341-w2)]. (105) 



Fig. 8. Normalized impulse responses of the linear system with 
its given cw squared-magnitude response in Eq. (90) and the 
approximate squared magnitude in Eq. (95). 



The procedures demonstrated in example 3 can 
be easily extended to cases with more resonant 
frequencies, where there will be more terms in Eq. 
(89). More algebraic processes will be involved. 
When extracting transfer functions from the 
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approximate squared magnitude, we will find more 
combinations for those with nonminimum phases 
while there is still only one solution for the transfer 
function with a minimum phase. Once the transfer 
functions are obtained, the remaining task for 
determining the corresponding phases and impulse 
responses is relatively straightforward. Determina- 
tion of passive elements and a specific circuit 
structure to represent the extracted transfer 
function is a typical network synthesis problem 
[4, 8, 9], which is not within our scope of analysis. 
Once a network is synthesized, we can then use this 
model to make further analysis and even measure- 
ment of the network response due to any excita- 
tion, cw or pulse, with little effort and cost. 

One caution, however, must be exercised for 
dealing with the cases where two resonant frequen- 
cies happen to be very close together. In these 
cases, the half-power frequencies associated with 
each resonant frequency must be entered into the 
computing process with smaller values than the 
actual values such that the final approximate 
squared magnitude still exhibits two distinct 
maxima at those resonant frequencies (rather than 
smeared together to have only one maximum). 

The three examples presented so far are simu- 
lations where the given magnitudes are expressed 
in terms of neat mathematical functions and the 
resonant frequencies are small and easily manipu- 
lable numbers. In the real world, this is definitely 
not the case. Our goal is still to deduce an approxi- 
mate squared magnitude from the given measured 
cw magnitude data so that a set of transfer 
functions and related characteristics can be deter- 
mined and analyzed. Another example under this 
situation will be presented later in Sec. 5. 



to use <uo = 2 to avoid manipulations with large 
numbers. After obtaining \H2a(J<^)\^> H^^is), 
0m(<w), and h'2a{t) in Eqs. (41) through (44), we 
can apply Eq. (108) with B =2^(10)^ to transform 
the results from <u to <u' with (o=o)'/B. Thus, the 
solutions in Eqs. (41-44) become respectively 



\G^ow'r = 



0.69315^ 



(o'*-8B^w'^+l6.6931B* 



G 0.8326g^ 

^^^ ^ s'^ + 0.414155 ' + 4.08575 ^ ' 

0™(<u') = tan-'[0.41415<u7('4.085752-<u'^)]. 



and 



(109) 



/i2,(0 = 0.414LB e-"-^'"sin(1.010750. 

These procedures apply also to the other type of 
second-order and higher-order transfer functions. 
With this explained, we now are ready to give 
another example based on the real-world data 
shown in Fig. 9. The data represent the measured 
but normalized electric fields (magnitude) of verti- 
cal polarization, reflected from a helicopter when it 
is irradiated by an impulse signal. By examining the 
curve in Fig. 9, we notice four significant resonant 
frequencies at 16.50, 26.25, 41.00, and 53.375 MHz. 
The frequency near 3 MHz is ignored because 
its magnitude response is rather small (close to 
background noise). It can, however, be added if 
necessary. 



5. Frequency Transformation 

In Sec. 4 the variable <u was loosely called 
frequency. Strictly speaking, <u is the normalized 
radian frequency. It can be translated into any 
frequency of interest by a simple frequency 
transformation [8], 



0)' = B<o , 



(108) 



where 5 is a normalization constant, <u is the 
normalized radian frequency, and <u' is the actual 
radian frequency. 

In presenting example 1 in Sec. 4.1, we cited 
(Uq = 2 as the resonant frequency. If the actual 
resonant frequency occurs at 20 MHz, we should 
have used (Uo = 4ir(10)^ rad. Instead, we chose then 
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Fig. 9. Measured electric-field magnitude (vertical polarization) 
reflected from a helicopter when radiated by an external 
impulse signal. 
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To simplify the analysis, we temporarily designate 
6>o 1 = 2ir X 16.50 = 103.6726, 6>o 2 = 52.5ir = 164.9336, 
6>o 3 = 82ir= 257.6106, and 6>04 = 106.75ir= 335.3650. 
From Fig. 9, we also see that the respective 
maximum responses at these, resonant frequencies 
are 14.00 (22.92 dB), 52.52 (34.41 dB), 13.00 (22.28 
dB), and 5.33 (14.53 dB). Their half-power 
frequencies are approximately 6>ii = 29ir = 91.1062, 
<ai2=36Tr= 113.0973; 6>2i=467r= 144.5132, 6>22=597r 
=1853540; 6>3i=767r=238.7610, 6>32=907r=2817433; 
and 6>4i = 104Tr= 326.7256, 6>42=108.67r=341.1770. 
Here the first subscript refers to the resonant 
frequency, the second subscript 1 refers to the half- 
power frequency on the left side of the respective 
resonant frequency, and the second subscript 2 
refers to the half-power frequency on the right side 
of the respective resonant frequency. The two half- 
power frequencies, 6>i2 = 36ir and fi>3i = 767r, are 
not real, but extrapolated for the analysis. Also, we 
later use B = 10*, in accordance with Eq. (108), as 
the transformation constant. 

Even though the condition eoi + eo2 = 2eoo is not 
exactly satisfied at either of the four resonant 
frequencies, we choose to use the type of second- 
order transfer function in Eq. (10a) for obtaining 
the approximate squared magnitudes. Using Eq. 
(23), we have for the first resonant frequency at a>o 1, 

yl 1 = 3.1435(10"), 61 = 1.0980(10"), ai = 21.5420, 

and 

|//i0'<»)p = 9.8816(10*)/[6>"- 2.1496(10") 6>2 

+ 1.2056(10^)]; (110) 

for the second resonant frequency at 6>o2, 

y42 = 3.5377(100, 62 = 2.8025(10"), 02 = 40.5358, 

and 

|//2(»P = 12.5153(10"')/[6>" - 5.4406(10")6>^ 

+ 7.8538(10^]; (111) 

for the third resonant frequency at 6>03, 

y43 = 1.4909(100, 63 = 6.7347(10"), 03=44.3546, 

and 

|i/3(;6>)|^ = 2.2228(10"')/[<i>''-l-3273(1006>' 

-1-4.5356(10'')]; (112) 



and for the last resonant frequency at (oo4, 

^4=2.5723(10"), 64 = 1.1257(100, 04 = 14.3871, 
and 

|i/40'<»)P=6.6166(100/[6>"-2.2494(1006>2 

+ 1.2673(lO"0],- (113) 

The final approximate squared magnitude is then 

\HUa>)\'=\H^ija>)\' + \H2(ja>)\'+\m(jo>)\' 
+ \H4ao>)\' 
= lA9Q3(W)N(eoyD (a>^) , (114) 
where 

iV(6>0 = 6>'2-3.6694(1006>"'+5.1348(10'0<i>'' 
- 3.4113(10'0 <^^ + lM2Q{lO^°)o}* 
- 1.3940(10^")6>2 + 6.2996(10") 

= [6>"-2.1747(10")6>2 + 1.2577(100] 

X [6>"-1.2098(1006>^ + 3.9741(10')] 

X [6>"-2.2422(100<»' + 1.2615(10'")] , 

(115) 

and D(6>^) is the product of the four denominators 
in Eqs. (110) through (113). The magnitude in 
Eq. (114) is shown in Fig. 10 together with the 
component magnitude functions obtained in Eqs. 
(110) through (113). Comparing Figs. 9 and 10, we 
see, except the frequency scale, the approximation 
in Eq. (114) is generally very good. The dominant 
features at 6>o2 and its half-power frequencies are 
indeed excellent. The shifts in 6>oi, and 6>o3, are 
minor. The position of 6>o4 remains practically the 
same. The only major changes are the magnitudes 
at 6> = and 6>o 1. This deficiency can be improved 
if we choose the second-order transfer function of 
Eq. (10b) or the approximate squared-magnitude 
function in Eq. (49) with c = at the beginning for 
|//20<")P- From Eq. (115) we already see the large 
coefficients even when we used the normalized 
frequency to begin with. If we wish to convert the 
frequency into megahertz, the numerator in Eq. 
(114) will become 

(l/5'O [<»''-2.1747(10")5^6>2 + 1.2577(10^5"] 
X [6>" - 1.2098(1005 ^6>^ + 3.9741(1005 "] 
x[6>"-2.2422(1005W-l-1.2615(10'05"], 

(116) 
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and the denominator of Eq. (114) will become 

(1/5'") [ft)*- 2.1496(10^5^6)' + 1.2056(10^)5 '] 
X [ft)* - 5.4406(10^)5 W+ 7.8538(10^)5 *] 
X [ft)* - 1.3272(10^5 W+ 4.5356(10')5 *] 
X [ft)*-2.2494(10^)52ft)2+ 1.2673(10"')5*] , 

where 5 =10*. 




20 30 40 

Frequency, MHi 



Fig. 10. The approximate magnitude to that in Fig. 9, together 
with component magnitude functions identified for each 
resonant frequency, where \H(Ja})\=[\Hi(_j'ai)\^+\H2U'^)\^+ 
|//,0•a.)|2 + |//40•<")l']"^Eq.(114). 

Referring to Eq. (114) and setting \H(j(o)\^ 
=H(s)H(-s)\s^ja, we obtain 

His)H(-s) = 1.4903(10")A^( -«')/£> (-«'), 

where ^ ^ 

N(-s^)=Ni(+)N2(+)N3(+)Ni(-)N2(-)N3(-), 

(119) 
and 

D(-s^)=Di(+)D2(+)D3(+)D,( + )Di(-) 

xD2(-)D,(-)D,i-), (120) 



with 



M( + )=s' + 26.1252s + 1.1215(10*), 
Ari(-)=s' -26.1252s + 1.1215(10*), 

N2( + )=s^ + 71A229s +6.3040(10*), 
N2i-)=s^-71A229s +6.3040(10*), 

Ar3( + )=«' + 20.5139s +1.1232(10^), 
A^3(-)=52-20.5139s +1.1232(V); 

Di( + )=s^ + 21.542Qs + 1.0980(10*), 
Di(-)=s2-21.542Qy + 1.0980(10*), 



D2( + )=s'+40.535&s +2.8025(10*), 
I>2(-)=s'-40.535&s +2.8025(10*), 

D3( + )=s^ + 44.3546s +6.7347(10*), 
£)3(-)=52-44.354fo +6.7347(10*), 

D^( + )=s^ + 14.3871s + 1.1257(10^, 
£)4( -) =s '- 14.3871s + 1.1257(10^). 

Since we require the system to be stable (no 
poles in the right half of the s -plane), we have to 
assign £)i( + )£)2( + )I>3(+)£>4( + ) as the denomi- 
nator for H(s). Thus, £)i(-)£)2(-)D3(-)I>4(-) 
belongs to H(—s). As far as the numerator for 
H(s) is concerned, we have many choices from 
Eq. (119). WhenM(+)A^2(+)A^3( + ) is assigned as 
the numerator of H(s), Ni(-)N2i-)N3(-) then 
belongs to H(—s). In this case, there are no zeros 
in the right half of the s -plane. The result is a 
minimum-phase transfer function. We then have 

H^is) = 3M04ilO=)Ni( + )N2( + )N3(+)/ 
[£).( + )£)2( + )D3(+)D4( + )] 

Fis+Gi F2S+G2 



. (FiS+Gi F2S+G2 

= 3.8604(100 -h + -^ — - + 

^l £),( + ) D2i + ) 

F3S+G3 F4S+G4 1 

£)3(+) "^ D,(+) r 



where 



(121) 



F, = 2.3245/10*, G, = 1.3664/10^ 
F2 = 5.3398/10*, G2 = 0.8742, 
F3= -6.8335/10*, G3 = 0.1055, 
F4=-8.3091/10^ G4 = 1.7752/10^ 

Applying the following two Laplace transform 
pairs [9] to Eq. (121): 

(s+a)/[(s + ay + fi^] <^ e-'-'cosfit, (122a) 
and 
l/[(s+a)'+^'] ^ (l//3)e-"'sin/3f , (122b) 

we obtain the impulse response of this system, 

MO = 3.8604(10'){(0.2325 cos /3,f +0.1071 sin /3,0 

X e-"" + (0.5340 cos )32f +5.1959 sin fiit) e'"^' 
+ (-0.6833 cos Pit +0.4665 sin ^^t) e""'' 

+ (-0.0831 cos )34r + 0.0071 sin fi^t) e""^' } , 

(123) 
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where 



ai = a ill, and /Si = V^P(«i^. ' = 1.2,3, and 4. 
More specifically, we have 



a, = 10.7710, 02=20.2679, 
04=7.1936, /3i = 104.2305, 
/33 = 258.5634, /34= 335.4422. 



03=22.1773, 
/32 = 166.1740, 



Equation (123) shows that at f = 0, hm{t) = Q. This 
agrees with the result predicted by the initial-value 
theorem [10, 12]. The largest coefficient is with 
sin^Szf associated with the second resonant 
frequency. This is obvious when we refer to Fig. 10 
where woz is dominant. 

With the impulse response so determined, the 
system's response to a general excitation can then 
be computed by convolution integral [12]. When 
referring to frequencies in megahertz, we simply 
modify the impulse response in Eq. (123) by multi- 
plying the coefficient 3.8604(10^), o,-, and 
/3,- (/ = 1, 2, 3, 4) by the transformation constant 
B = 10*. The impulse response hm(,t) before apply- 
ing the frequency transformation is presented in 
Fig. 11, where the coefficient of 3.8604(10^) in Eq. 
(123) has been dropped. We see, from Fig. 11, that 
the period is about 0.019 s, giving 0.019 /Sz^tt. The 
major maximum occurs approximately at 
fi = 0.008 5 with ;tm(f,)/3.8604(10=') = 5.3094, and 
the second maximum occurs approximately at 
^2 = 0.026 5 with A„(f2)/3.8604(10=')= -3.2652. The 
ratio of Am(^i)/|Am(^2)| = 1.6261, which is close to 
ga2((2-'i)_ Thus, the second resonant frequency is, in- 
deed, the dominant one [14]. 






0.04 
Time, s 



Fig. 11. Normalized impulse response of the linear system 
whose approximate magnitude of the transfer function is shown 
in Fig. 10. This is for the minimum-phase case with its transfer 
function given in Eq. (121). 



The associated phase is given by 

9m{0)) = 0l(ft))-|-02(<w) + 03(<w)+M<y) 

- 05(ft)) - 06(<w) - eii(o) , (124) 

where the first four component phases are due to 
I>i( + )I>2( + )I>3( + )I>4( + ), and the last three are 
due toM( + )A^2( + )A^3(+). That is, 

, r 21.5420ft) 1 

0i(co) = tan"' , 

^ ' L 1.0980(10") -ft)^ J 

, r 40.5358ft) 1 
W = tan-[^^^^^^^,^_^,J, 

, r 44.3546ft) 1 

62(0}) = tan"' : , 

^ ^ L 6.7347(10") -ft)^ J 



, r 14.38716) 1 

4(0) = tan"M z 

^ ' L 1.1257(10^) -6)2 J 



1257(10^^ 
26.12526) 



. r 26.12526) 1 
^-^") = ^^" L.1215(10")-ft)- J' 

, r 71.42296) I 
M-) = tan-[ ,J, 



and 



07 (ft)) = tan' 



. 6.3040(10") -ft)^ 
20.51396) 



, r 20.51396) I 
L 1.1232(10^) -ft)2 J' ^^^^ 



Expressing the phase in terms of megahertz, we 
multiply the numerator inside the arctangents by 
the normalization constant B and the constant term 
in the denominator by B ^. The minimum phase 6m 
in Eq. (124) before frequency transformation is 
presented in Fig. 12. 

Seven other possible solutions for the transfer 
function with nonminimum phases can be obtained 
from Eqs. (118) and (119) as 

H„i(s) = CNti + )N2( + )N3i-), 

Hnlis) = CNl( + )N2(-)N3( + ), 

H,^(s) = CNi(-)N2( + )N3( + ), 
HnA(s) = CiVi( + )A^2(-)iV3(-), 
Hr^is) = CN,i-)N2i + )N,i-), 
H„e(s) = CNi(-)N2(-)N3( + ), 

and 

H.i(s ) = CNi(-)N2(-)N3i-), (126) 

where 

C = 3.8604(10y[Dii + )D2i + )D3i + )D,( + )]. 
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Fig. 12. Phases of the sample linear system whose approximate 
cw magnitude response is shown in Fig. 10. 



Following the same procedures of partial 
fractions as in Eq. (121), we list the expansion 
coefficients for the 7 nonminimum-phase cases in 
Table 3. 

The corresponding impulse responses also take 
the same form as in Eq. (123) with the same a,- and 
Pi but with different coefficients associated with 
the cosine and sine terms. These coefficients are 
listed in Table 4. 



These impulse responses with nonminimum 
phases, without including the frequency transfor- 
mation, are shown in Figs. 13 (a) and (b) to com- 
pare with that in Fig. 11 for the minimum-phase 
case. It happens that the first maxima ofhai(t) are 
all below the first maximum of hm(t), and that the 
first nulls ofha\(t) are also closer to the origin than 
the first null of /im(0- More energy is concentrated 
near t = in hm(t) than any of the hm(t), 
i = l,2 7. 
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Fig. 13 (a). Normalized impulse response of the helicopter with 
the approximate magnitude response given in Fig. 10, but with 
the nonminimum-phase transfer functions given in Eq, (126). 



Table 3. Partial fraction expansion coefficients for nonminimum-phase transfer functions 



Non-Min. 


f 1 X 10' 


GiXlO 


F^y.Vfi 


G^ 


fjXlO' 


G3XIO 


f 4X 10^ 


G4XIO 


nl 


0,2290 


0,1470 


0,1253 


0.8842 


-0.8235 


0.6849 


4.6925 


0.1281 


n2 


0,2017 


0,2034 


-3,0480 


0.8467 


2.8954 


1.1708 


-0.4910 


0.2434 


n3 


-2,5895 


-0.0188 


3.2415 


0.8940 


-0.5711 


1.5046 


-0.8103 


0,0653 


n4 


0,1952 


0.2124 


-3,5248 


0.8083 


2.9958 


2,9139 


3.3377 


-1,2197 


n5 


-2.6114 


-0.1339 


2.8849 


0,9406 


-0.7427 


1.2150 


4.6927 


-0.1422 


n6 


-2.6634 


-0.8067 


-0,2450 


1.1784 


2.9444 


-0,6197 


-0.3600 


0,2700 


n7 


-2,6541 


-0.9251 


-0.8178 


1.1789 


3.2050 


1.0964 


2.6693 


-1.4034 



Table 4. Coefficients associated with cosine and sine terms in impulse responses for nonminimum-phase cases 



Non-Min. 


cosiSi/ 


sin/3ir 


co&^zt 


sinj82' 


cosftr 


sinftr 


cos;34' 


sin/34< 


nl 


0.2290 


0,1174 


0.1253 


5.3059 


-0.8235 
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Fig. 13 (b). Normalized impulse response of the helicopter with 
the approximate magnitude resfwnse given in Fig. 10, but with 
the nonminimum-phase transfer functions given in Eq. (126). 
Continuation of Fig. 13 (a). 

The phases associated with the transfer function 
in Eq. (126) are respectively: 

0nl = 01+02+03+04-05-06+07, 

0n2 = 01+02+03+04-05+06-07, 

0„3 = 01+02+03+04+05-06-07, 

0„4 = 01+02+03+04-05+06+07, 

0n5 = 01+02+03+04+05-06+07, 

0n6 = 01+02+03+04+05+06-07, 

and 

007=01+02+03+04+05+06+07, (127) 

where 

0,, 1=1,2, ..., 7, are given in Eq. (125) 
before normalization. 

Graphs for 0ni are also plotted in Fig. 12 for 
comparison purpose. Clearly we see that 0ni > 0m, 
because each component phase given in Eq. (125) 
is nonnegative, varying from to tt as <u varies from 
to 00. 

6. Consideration of Energy Contents 

To assess the ability of a system to withstand 
damage from an external unwanted excitation, it is 
often useful to compute the energy content associ- 
ated with an impulse response [15]. Indeed, if h(t) 
represents a voltage waveform across a 1 fi resis- 
tor, the quantity 



E = f h\t)dt, 
Jo 



(128) 



equals the total energy delivered to the resistor by 
the impulsive excitation [2]. Equation (128) also 
represents the area under the cuTwe h^(t). 

The energy E may also be computed, in view of 
Parseval's theorem [2], by 



E=^f_jHU<ord<o. 



(129) 



Thus, when the minimum-phase impulse response 
/ira(0 Slid the associated nonminimum-phase 
impulse response hn(t) have an identical \H(j'io)\^, 
their respective total energies [in 0^/=S <»] are 
equal even though hn(t)<hm(t) during the initial 
period near / = + , as discussed in Sec. 4. These 
facts can also be demonstrated by referring to the 
examples given earlier. 

For example 1 in Eq. (40), \H2aU^)\^ '^ given 
in Eq. (41). Its energy content, according to 
Eq. (129), is 



0.6931 r 
Itt J- 



d(o 



(0^-8(0^ + 16.6931 



0.6931 r _ 

IT Jo ( 



d(o 



(o^+p^X(o^+p*^) 



0.6931 f° r_J I ] 

~;1.66507r Jo l(o^+p^ ~ (o^+p*^ \ 



= 0.2049, 



(130) 



where 



p^ = -4-; 0.8325, and 

p* = complex conjugate of p. 

The last step in Eq. (130) is obtained by using 

cos(qx)dx TT 



1 



x^+p^ 2p 

with q^O 
and Re(j>)>0. 



= — Q-tP, 



(131) 
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We can also obtain the energy by referring to 
the impulse response given in Eq. (44) for the 
minimum-phase case, 



£»= f hl(t)dt 
Jo 

= 0.1715 f e-'''"*'sin2(2.0107Odf 
Jo 



= 0.2049, 



(132) 



which is indeed the same as in Eq. (130). 

Using the corresponding impulse response given 
in Eq. (47) for the nonminimum-phase case yields 
the same result. That is. 



JO 



h'n{t) dt = 0.2049. 



(133) 



For example 2 presented in Eq. (75), the corre- 
sponding \H2b 0<^)|^ ^m(0> ^nd h„{t) can be found 
respectively in Eqs. (76), (81), and (82). The total 
energy for this system is 



E=:^ j_. \f^^»Uo^Wdc 



-f 

Jq 

- f 

Jo 



hUf) dt 



hl(f) dt 



= 2.3124. 



(134) 



If we replace the upper integration limit « in 
Eq. (128) by a finite T, we can analyze the energy 
content absorbed by the system during the initial 
period after an external excitation is applied. 
Referring again to example 1 with /i2a (0 given in 
Eq. (44) and hn(t) in Eq. (47), and carrying out the 
details, we have 

hl(t) = 0.1714 e-""sin2^^ 

= 0.0857 e-""(l- cos 2/30, 



h^(t) = 0.1270 e-2' - (0.2540 cos pt 
+ 0.1950 sin /gOe-'-'""*" 
+ (0.1009 + 0.0261 cos 2/3f 
+ 0.0975 sin 2/306""", 

£«,= f /iL(Odf =0.2048 -0.2070 e""'^ 
Jo 

- 0.0052 e-'"^( - ai cos 2)37 
+ 2/3 sin 2)37), 



and 



£«= I /i„'(Odf = 0.1413 + 0.0635(1 -e-2^) 
Jo 

+ ( - 0.2437 - 0.0246 cos 2/37 

+ 0.0039 sin 2/37) e""'^ 

+ (0.1270 cos /37 

-0.0501 sin /37) e-'-^"™^, (135) 



where 



cii = 0.4141, and /3 = 2.0107 



Numerical results for both Em and E„ are shown in 
Fig. 14 and indicate clearly that £„<£„ [12]. They 
are equal only when 7 -» oo . This reconfirms that 
the impulse response and transfer function with a 
minimum phase deduced from a given magnitude 
can be used as the worst case for analysis purpose, 
as far as the initial impact to the system under 
study by an external unwanted source is concerned. 






10 
Time, s 



Fig. 14. Energy contents of the sample system whose transfer 
functions are given in Eqs. (42) and (45). 
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For the practical example shown in Figs. 9 and 
10, we present the results on energy content in 
Fig. 15. Again, we have Em <Em , / = 1, 2, . . . , 7. 
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Fig. 15. Energy contents of the helicopter whose approximate 
cw magnitude response is given in Fig. 10. 



7. Conclusions 

We have used a simple method known in classi- 
cal network theory to determine the complete 
characteristics for an unknown linear system from 
a given cw magnitude response only. These char- 
acteristics include possible different transfer 
functions, their phases, and the corresponding 
impulse responses. Only one transfer function is 
minimum phase. The main achievement is to 
deduce an approximate squared-magnitude func- 
tion in the form of a ratio of two even polynomials 
based on the outstanding features in the given 
magnitude response, such as resonant frequencies 
and bandwidths. The remaining procedures for 
obtaining the complete system characteristics are 
exact. Four examples, three simulations and one 
using measured data, have been given to illustrate 
the proposed method. We have written software 
that greatly facilitates application of this technique. 
It first performs numerical calculations necessary 
to obtain the system transfer functions from the 
measured magnitude-frequency input data, and 
then gives impulse responses. We also have shown 
that the minimum-phase case, through its asso- 
ciated impulse response and energy content, 
constitutes the most pessimistic estimate as far as 
the initial threat to the system is concerned. 
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